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The one dimensional quantum walk of anyonic systems is presented. The anyonic walker performs braiding 
operations with stationary anyons of the same type ordered canonically on the line of the walk. Abelian as 
well as non- Abelian anyons are studied and it is shown that they have very different properties. Abelian anyonic 
walks demonstrate the expected quadratic quantum speedup. Non-Abelian anyonic walks are much more subtle. 
The exponential increase of the system's Hilbert space and the particular statistical evolution of non-Abelian 
anyons give a variety of new behaviors. The position distribution of the walker is related to Jones polynomials, 
topological invariants of the links created by the anyonic world-lines during the walk. Several examples such as 
the SU(2)jt and the quantum double models are considered that provide insight to the rich diffusion properties 
of anyons. 

PACS numbers: 05.30.Pr, 05.40.Fb 



I. INTRODUCTION 



Quantum versions of random walks came to prominence in quantum information theory through the search for 
efficient new quantum algorithms. Since classical random walks provide the techniques for some of the best classical 
algorithms, it was natural to look for a quantum equivalent. Typically, a quantum walk gives a quadratic algorithmic 
speed up over classical ones JU and for some problems an exponential speed up is possible [2]. This is a generalization 
of the faster spreading behavior shown by the simplest quantum walk, a single walker on an infinite line. 

Quantum walks have since found wider applications than their algorithmic origins. In essence a quantum walk is a 
discrimination of a diffusion process. For example, they has been employed in modeling transport of charge or energy 
in biological systems [3], in physical systems to show derealization and quantum coherence 0,H,0], as a model 
system interpolating between quantum and classical behavior 17D, as well as a model amenable to quantum simulation 
of its asymptotics [80 . On the algorithmic side they have been proved to be universal for quantum computation [9, 10] 
and they have been physically implemented, both as a computation on a quantum computer [11] and a physical walk 
(e.g. lll2l[l3lD . Quantum walks have been studied theoretically on many different structures: in higher dimensions 1 14], 
under the effects of decoherence lfl5ll . with Dirac spin particles [16], and with a huge r ange of variations on the basic 
walk dynamics. For accessible introductions and reviews, see, for example, lfl7i [18111911 . 

Here we focus on the question of how the quantum walk distribution can be affected by the statistical properties of 
the walker. Non-interacting bosons in quantum walks have effectivelyalready been shown to have identical properties 
to single quantum walkers through experiments with coherent light lr7L [1211 . Omar et al. lf20ll studied quantum walks 
with two bosonic or fermionic walkers. Entanglement in the initial state of these walkers gives a strikingly different 
behavior with respect to the standard quantum walk. 

More exotic statistics than the bosonic or the fermionic can be found in two-dimensional systems where anyons can 
appear. Exchanging two anyons can introduce phase factors or unitary transformations. These evolutions are different 
representations of the braiding group. Envisioning how quantum walks can be performed with anyons one quickly 
realizes that the problem takes the form of diffusion of a walker-anyon in the presence of other regularly arranged 
static anyons. In particular, we embed a one dimensional quantum walker into a two dimensional medium with a 
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canonical line ordering of the anyons from left to right. When the anyonic walker is moving among the static anyons 
then anyonic braiding is realized evolving the overall state of the system in a non-trivial way. Hence, by considering 
an anyonic walk we enrich the one particle diffusion problem with statistical properties. 

Anyonic quantum walks is not only an academic curiosity, but it can be considered as a quantum simulator modeling 
the diffusion of extended objects. In order to keep track of the statistical properties of anyons it is useful to visualize 
their world-lines that encode their braiding history. When the walker moves from its initial position it braids with the 
straight world-lines of the static anyons. Finding the walker at a certain position after a number of steps requires the 
consideration of all the different paths of the same length with the same initial and final positions. This gives rise to 
world-line links the complexity of which increases fast with respect to the number of the walker steps. The diffusion 
of extended objects has applications in biophysics and polymer physics 02111 . Here, we are interested in the quantum 
version of such diffusion processes. 

Before turning to the anyonic quantum walks we briefly introduce the standard quantum walk as well as the anyons 
and their properties. In Sec. [II] we present in detail the general formalism of the anyonic quantum walk. In SecHIIlwe 
demonstrate that Abelian anyonic quantum walks have similar asymptotic behavior as the standard quantum walks. 
In Sec.[lV]we consider the non- Abelian anyonic quantum walks and we express the walker probability distribution 
in terms of the Kauffman brackets. Several examples are explicitly solved for a small number of walker steps that 
indicate a wide range of diffusion properties ranging between the classical and the quantum walks. Finally, in SeclVl 
we present our conclusions and we point out interesting future directions. 



A. Quantum walk on a line 



Quantum walks have been studied in both the continuous time l22ll and discrete time ll23l l24ll versions. We briefly 
describe the discrete time quantum walk on the infinite line, which will be later generalized to the anyonic case. It is 
defined in direct analogy with a classical random walk: there is a walker carrying a coin which is tossed each time step 
and the walker moves left or right according to the heads or tails outcome of the coin toss. 

We denote the basis states for the quantum walk as an ordered pair of labels in a "ket" |s) spa ce <8> |j)spin> where seZ 
is the position and j S {0, 1 } is the spin-like state of the coin. A unitary coin operator is used at each time step and then 
a shift operation is applied to move the walker to its new positions. The simplest coin toss is the Hadamard operator 
H, defined by its action on the basis states \s, j) as 

H\s,0) =^=(|*,0>+K1» 

ff|*,l) =^(|*,0>-kl», (i) 
and the shift operation T acts on the basis states as 

T\s,0) =\s-hQ) 

T\s,l) =\s+l,0). (2) 

A single step of the quantum walk consists of W = TH applied to the quantum walker plus spin. After performing t 
steps of the quantum walk, where the walker is initially in state | ^(O)), we obtain the final state 

|*(0) =W r I O(0)). (3) 

Ultimately we are interested in the probability distribution P(s,t) of the spatial location, s, of the walker at time t that 
is given by 

(P(l, t),P(2,t), ...) = diag (tr^Q^f)) <<*>(*) I] ) . 

As the walk progresses, quantum interference occurs whenever there is more than one possible path of t steps that 
leads to the same position. This interference is both constructive and destructive, which causes some probabilities to 
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be amplified or decreased at each timestep. This leads to the different behavior compared to its classical counterpart. 
In the latter the position of a walker, following a classical random walk on a line, spreads out in a binomial distribution 
about its starting point. In Refs. 12411 . Ambainis et al. proved that the quantum walk on a line spreads in 0(t 2 ) 
compared to a classical random walk which spreads in 0(t), where the spreading is being measured as the second 
moment about a single initial starting point. 



B. Anyons and their properties 



It is commonly accepted that point-like particles, elementary or not, come in two species, bosons or fermions. This 
statistical label is determined by the behavior of their wave function when two identical non-interacting particles are 
exchanged. For bosons the wave function is unchanged while for fermions it acquires a phase of Jt, i.e., changes 
sign. These are the only observed statistical behaviors for particles that exist in our three dimensional world. If one 
restricts to two dimensions there are more possibilities in the particle statistics. In this case, when two particles are 
exchanged their wave function can be evolved by an arbitrary phase factor or even a unitary operator that creates 
superpositions in an internal space of the particles. These particles, named anyons by Frank Wilczek, can appear 
as effective quasiparticles, such as vortices, in two dimensional many body systems. They are manifested in the 
Fractional Quantum Hall Effect (FQHE) 125[ 12611 and they are expected to ap pea r in p-wave superconductors B27L 12811 . 
topological insulators HQ and a variety of lattice models JllIIlSBHl]. 



Anyons can be created from the vacuum in pairs, they can be braided around each other and they can be fused 
together. While the pair creation results in a particle and its antiparticle much in the same way as usual particles, 
braiding operations and fusion of anyons are rather unique. This is due to the non-local degrees of freedom carried 
between anyons that give rise to their exotic statistics. One can think of it like an internal 'spin'. When two anyons 
are braided then their internal 'spin' state either acquires a phase factor, e li> (Abelian anyons) or it gets evolved by a 
unitary matrix, b (non- Abelian anyons). The fusion corresponds to the tensor product of these 'spins' with a restriction 
on the maximum possible 'spin' value. When two anyons, a and b, are fused then more than one outcomes are possible 
depending on their internal state. Symbolically we have a x b = c + d + .... 




FIG. 1: The (a) F and (b) R moves dictated from the fusion and braiding properties of anyons. 



Three anyons a, b and c can be fused to a final anyon d in two different ways. Either one can fuse a with b and its 
result i with c or one can fuse b with c and its result j with a. These two different processes correspond to two different 
bases of the collective state of the anyons. The different ways of combining particles a,b,c to yield d is given by the 
recoupling formula: \{ab)c — ► d;x) = EyC^Ac)* l( a (^ c ) ~~ * d;x') as seen in Fig. 02 a). The braiding of two anyons a, b 
with fusion channel c is described by R c ab as seen in Fig.QJb). To obtain the braiding result of two anyons that do not 
have a direct fusion channel (see for example anyons b and c on the left hand side of Fig. [TJ a)) one can use the F matrix 
to obtain their braiding element, F~ l RF. Anyons also carry a quantum dimension d. This accounts for the dimension 
of the internal space of the anyons and provides the scaling of the total Hilbert space. For example, n anyons with 
quantum dimension d have a Hilbert space with dimension d". Abelian anyons have d = 1 and non- Abelian ones have 
necessarily d > 1. Note that, in contrast to usual spin, d can be an irrational number. For more details on anyons and 
on topological quantum computation see, e.g. 13611 and references therein. 
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II. WALKING ANYONS 



In order to reveal the statistics of anyons in a quantum walk, we consider a single anyonic walker braiding around 
others in fixed positions. The anyons involved are all of the same type a. The walker has an attached die degree of 
freedom which gives it distinguishability. However this distinguishability does not negate the action of braiding on 
the non-local statistical degrees of the system. The anyons are placed in a canonical order in the plane with position 
labelled by an integer s, as seen in Fig.|2fa). For an anyon at position s, other anyons with positions s 1 < s (V > s) will 
be said to be on the left (right). We restrict to n mod 4 = 2xodd in order to have an odd number of pairs such that the 
initial state with the marked anyonic walker and its partner located at the middle will have an equal number of pairs to 
their right and left. This facilitates the construction, if desired, of a symmetric evolution. 
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FIG. 2: The anyonic quantum walk, (a) A quantum walk of the anyonic walker with a four state die attached. It walks along one 
of four braids depending on the state of the die on sites labelled by s = 1 to s = n (white circles) positioned in between stationary 
anyons (blue circles), (b) A basis state \^if(a\ ,d2, ■ ■ ■ ,a«-2)) of the non-local fusion space representing the sequential topological 
charge measurements of n type a anyons with total charge zero (vacuum, also denoted as 1). 



Basis states are labelled by ^V{a\,a2, ■ ■ -fln-2)) according to the fusion outcomes as indicated by Fig. |2|b). The 
state with nearest neighbor vacuum pairs depicted in Fig. [2J a ) is 

where 1 is the vacuum particle, as can be verified by considering the consistent fusion ordering. If desired one could 
modify the protocol to give a symmetric distribution for the walker by changing the initial state to provide for a 
symmetric spatial distribution for the anyons. Taking the vacuum partner of the walker and winding counterclockwise 
halfway around all its neighbors, the initial fusion state l^o) = b n b n -\ ■ ■ ■fr( n +3)/2l v ^'(lj cr i 1 , CF, - - - , 1, cy)) provides for 
the walker anyon to have the same number of vacuum pairs on its left and right. 

The Hilbert space of our system is 

# («) = Canyons («) ® #die(«)- (4) 

The die contains both a discrete spatial index indicating its position at the location of the anyons on the plane and 
4 internal spin states: #die(") = ^space («) ® ^4pin, where # spaC e = s P an c{l s )}"=i and ^spin = s P an c{|;)}]=o- The 
anyonic Hilbert space is further decomposed as Canyons («) = ^iocal(w) <8 ^fusion («) the first factor describing degrees 
of freedom local to each anyon and the later factor the non-local or fusion degrees of freedom. We restrict to a 
contractible surface so that the homological degrees of freedom are trivial. Furthermore, the local degrees of freedom, 
while measurable, will be conserved during this walk and will be ignored. 

Fusion degrees of freedom are determined by the fusion rules a,- x aj — Y^kNafajak, where N c ab £ N counts 
the number of ways to combine anyons of type a and b to obtain c. The corresponding Hilbert state space 
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of n o anyons with total charge zero is #f U sion(") = span c {| x P({a_ / -}-"~j ))} with dimension dim(.?/f us j on («)) = 
£ B] „ 2 0n _ 2 NaaNai a Nal a ■ ■ ■ Nl a . We will restrict to cases where the fusion spaces are one dimensional, i.e. N£ b < 2. 
Braiding relations are determined by the F and R matrices. The action of the generators of the braid group 2J„ 

on ^fusion (n) is given by the amplitudes 

'n\',,«,Xa, k=l 

m=k 

ff^^IC^a]- 1 )^ 1 ^^^,. k = 2 

m=k x 

ff S 4* II S^La^aar 1 )?" Xa(^ 2 aa)^, 2 < k <n- 1 ' 

f=l m=A: x 

n-2 

n § <,«Xa 3 . *=«-l 

where the formula b = F~ l RF is employed when the walker does not have direct fusion channels with the rest of the 
anyons. 



(V({a'})\b k \V({aj})) 



The initial state of our system is 



|«&(0)) = l^fusionh = ^ ) space | ¥> spin- 



The quantum walk algorithm consists of an iterative sequence of operations W — TU where U is the die tossing 
operation that acts trivially on the position of the walker and T is the conditional braiding operator. We can introduce 
a virtual tensor product structure to the spin states: |0) sp i n = |0) x |0) y , 1 1 ) sp in = |0) x |l) y , |2) sp i n = |l) x |0) y , |3) sp j n = 
|l) x |l)y where |0) y (|l) y ) are modes which are behind (in front of) their neighbors, and |0) x (|l)x) are modes moving 
left or right. For the die toss we pick (in the product basis) 

which can be interpreted as a pair of beam splitters on the two axes. Each component of U is a symmetrized version 
of the Hadamard operation of Eqn. ((TJ. Of course a continuum of tossing operations could be considered but this 
decomposition into a virtual tensor product has advantages as described below. The conditional braiding operation T 
has the following action 

7 , | X P)fusionk)space|K) x |Y)y = ^k_T V) fusion I ^ + 2K - l)space|K) x |Y)y 

over the spatial range < s < n + 1 . 

It is convenient to introduce the unnormalized spinor 

\$(s,t)) = { spm {0\^(s,t)),^ m {l\^( S ,t)), sp - m (2|4>( i ,f)), spin (3|<I>(.,f))) r , 
so that an arbitrary state at time t is 

l*C)> = 1 !*&*)>■ 

In order to avoid boundary conditions on the walk we study the behavior for t < n/2. 



III. ABELIAN ANYONIC QUANTUM WALK 
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A. Variance of Abelian anyonic walks 



Employing the method of Brun et al. IL37L l38i l39ll one can analytically evaluate the variance of the distributions 
corresponding to the Abelian anyonic walks by determining the moments of the resulting distribution. For simplicity 
we relabel the positions so that the initial walker position in the middle of the anyonic chain is given by sq — 0, i.e. 
|4>(0)) = \s = 0) sp ace <£> | V) spin- In Fourier components we have 

"=<»=/-,>■ 

In the Fourier transformed basis the evolution operator W acts as W(\k) (8 |\|/)) = \k) ®M k \\\i) so after t steps the state 
evolves to 



One can write a closed form for the moments of the resulting distribution. By employing the expansion of the 
identity in the position s basis, l = £ s \s) (s\ and (s\k) = e~' ks we obtain 

= ^l5'Y^/^' S( "^*(°)lK)'(^')'l*(0)> 
;m r r 

= ^J dk J dk'8^(k-k')(<P(0)\(Ml)'(M kl r\<P(0)). (5) 
Now we can integrate by parts the derivatives that act on the 8-function. For the case of m = 1 we obtain 

We have that dM k j dk = —i(Pg — Pi)M k , where Pg^ are the projector operators for moving the walker right or left with 
Pr+Pl= 1 ■ Finally, we obtain 

= ^ t / dk(mwt)ifa-PL)m J \*m- 

Similarly for ra = 2 we have 

(A = i 1 1 I dk(<P(o)\(Mly(p R -p L )(M k y-j\p R -p L )(M k y)<P(o))- P) 
zn j=lj'=\ J 

One can further employ the eigenstates, \Xi(k)), and eigenvalues, X/(k), of M k (I runs through the eigenstates) to 
explicitly evaluate the moments. Let us expand the initial state in this basis 10(0)) = T,i c i{k)\^-i(k)). After t steps the 
state becomes (Ma)'|4>(0)) = Y.I e ' X '^'ci(,k)\Xi (k)). With this in mind the first moment becomes 



(s) t =t-l f d*£c?(*M*)<M*)lftlM*)> £ ^ 

" J i n ;— 1 



J[X,,(k)-X,(k)]j 

U' j=i 



If the matrix M k is non-degenerate then most of the terms in (s) t will be oscillatory and they will average to zero over 
time. The non-zero contributions comes from the diagonal terms 

(s), = {l--J dkY^\ci{k)\ 2 {h{k)\P L \h{k))y + 0&ci\\?Ao^ terms. (8) 



Similarly for the second moment we have 



(A= l- 



\ f dk H \ c i( k )\ 2 ( k i ( k ) \ p L\h i k )) (h ( k ) \ p R\^i( k )))t 2 + Oscillatory terms and lower t orders. (9) 



For the Abelian anyonic walker, 

M k = e ik M++e- ik M- = e -i^x^ZxmZy e il (x x +X y )^ 
The eigenvalues of Mf, come in conjugate pairs 

Xi(jfc) = e- ! 'P-«, X 2 (k) = e' M<0 h(k) = e-^ k \ X A {k) = e^ (k \ 



where p±(£) = cos 



-l 



\ ( cos(fc) cos(0) ± ^(cos 2 ((j))-2)cos 2 (/fc)+2 



(10) 
(11) 



The eigenvectors are 



|Xi(Jt)) = |v(P_)), \X 2 (k)) = -iY®Z\X 1 (k)), 



Mk)) = |v(P+)>, \X 4 (k)) = -iY®Z\X 3 (k)), 



where 



|v(P)> 



(g -2/(*+P) + ! + g -2iP _ g i(-fc+*-P) + e 2,(0- P ) _ 2e «(-t+*-3P) _ e ^+*-P))|0),|0) v 

+i( e -i«'+*+P)- e -2^+P) _ i + < r i (*-*+P))|o),|i), 

+,■(£-''(*+*+?)_ 1 _ e -2ip + e i(fc+t>-p))| 1 ) jc | ) y + ( e 2/((|.-P)_ 1)|1) X |1)J , 



(12) 



(13) 



with 1 / V^Y the normalization factor. Employing the eigenvectors and eigenvalues of Mk and using Eqs. (O, © the 
variance, v = (s 2 ) — (s) 2 , can be evaluated, as seen in Fig. [3] 
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FIG. 3: The variance, v = (j 2 ) — (s) 2 , of the Abelian quantum walk as a function of the number of steps, t, and the statistical angle, 
<]). The two plots shown are the analytical results (bottom surface), based on Eqn. {8} and {5), and the numerical ones (top surface) 
giving almost identical results. A quadratic behavior of the variance is witnessed as expected from the standard quantum case. 



B. Results 



It is rather straightforward to perform a numerical simulation of the Abelian quantum walk on the line for relatively 
large number of steps. This gives the exact distribution of the walker and one can obtain the corresponding variance, 
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v = (s 2 ) — (s) 2 , as a function of the number of steps, t, and the statistical angle, 0, of the anyons. The numerically 
obtained variance and the analytical results of the previous subsection are depicted in Fig. [3] appearing to be almost 
identical, with a small difference due to the neglect of terms in the analytic calculations. We observe the expected 
quadratic behavior that corresponds to the standard quantum case. There is a weak dependence on the anyonic statis- 
tical angle 0. This is due to the form of the coin evolution that includes the term e li>z >® Zx (see Eq. ([Toll). For the case 
where is a multiple of K /2 this term has a tensor product structure that splits the system into the tensor product of two 
quantum walks, each one with a two dimensional coin. Thus, we expect to obtain the same variance as the standard 
quantum walk with a single two dimensional coin. Note that these cases include the quantum walks of bosons (0 = 0) 
and fermions (0 = n). For generic values of the walk entangles the two coins. As a result the interference effect 
that gives the speed up of the quantum evolution gets diluted over the four dimensional space of the two coins and the 
slope of the variance is decreased. This is in agreement with previous studies, where quantum walks with multiple 
coins were considered ll37l l38ll . 



IV. NON-ABELIAN ANYONIC QUANTUM WALK 



In this section we would like to perform a quantum walk with non-Abelian anyons. For simplicity we slightly change 
the quantum walk. Now the walker moves to the left and to the right by always braiding in the anti-clockwise sense 
with its neighboring anyons. This can be achieved by employing only one coin (which has two degrees of freedom). 
One can show that for a given non-Abelian anyonic model the braidwords that result from this quantum walk have the 
same computational power as arbitrary braidwords. Initially, note that having static n anyons and braiding one moving 
anyon clockwise or anticlockwise around them can realize the full computational power of the anyonic model [40]. 
Restricting only to anti-clockwise braidings does not restrict the possible obtained unitaries. Indeed, for the finite 
braiding groups we consider here one can show that there exists an integer k such that for a generic braid element bi 
it is b\ = 1 => b^ 1 = bj l (see l4lll ). Thus, braidings equivalent to clockwise moves can be produces given enough 
anti-clockwise ones. Notably, fixing the chirality of the braiding means that Abelian sectors of any physical theory 
will not contribute to the walker's distribution since all walking paths accumulate the same overall phase. Hence the 
single coin protocol probes the truly non-Abelian statistics. 

Consider t steps of the quantum walker and, for brevity, define 

(-3 

B 'a = Y\ b s,+2(a l + ...a,- r )-(t-r)bs,+2a l +a 2 -2b s ,+a l -\- 
r=0 

The evolution of the system can be given by the unitary operator 

r 

V '= E E Y\Pa r U®B , a ®\s t +2{a l + ...+a t )-t){s t \. (14) 

s,e[H-l,n-f+l]3e{0, !}/ = ! 



Starting from the walker in position m and the coin qubit in state one obtains after t steps the state 

|¥(0>=V»|oc)|m> =E(fl^)lv)4l«)l4): d5) 

a r=\ 

where a is a t dimensional vector with r-th component or 1 depending if the walker is moving, respectively, left 
or right at the r-th step, h' s — m + 2(a \ + ... + a t ) — t and the initial fusion state of the anyons is |a) = |0101 ...01) as 
in 1431 . One can write the corresponding density matrix p' = | l P(f))( l P(f)| as 

P(0 - E(ri^^)lv)(¥l(Il^^)4l«)<a|4' t l« 3 )(4l- (16) 

a,a' r=l r=l 
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Importantly, the coin, walker position and anyonic Hilbert spaces can be traced out independently. For brevity let us 
set 

r=\ r=l 

and 

^=4I«)(«|4' + - 

We would like to obtain the position distribution of the walker regardless of the coin and anyonic states. For that we 
need to trace the coin states and the fusion space of the anyons resulting in 

trcointranyonp(f) = £ trej^try^ %) (n r s , I . 

Due to the choice of coin manipulations one needs to have paths a and a' that end up at the same coin state for the 
trace to be non-trivial. To simplify the presentation we consider the composite loop aa' corresponding to paths a and 
a' that start at the same point and end at the same point after the same number of steps. Thus, they are contributing to 
the diagonal part of tr coin tr anyon p(f ). 

For a Hadamard coin flip operation U =H (see Eq. (Q]i) one can easily evaluate that 

tr*4r = ^(-l) Z : 

where z is the number of all successive pairs of l's in the loop aa' . Hence, successive zeros and ones in the loop aa' 
can be compressed. While the coin trace is rather straightforward care needs to be taken when one takes the trace over 
the anyonic degrees of freedom. It can be shown [42] that the trace of unitary representations of the braiding group 
corresponds to the Markov trace of the braid graphs described above. As initial condition we adopt the anyons to be 
created pairwise from the vacuum. Hence, the trace of the braidings gives the Kauffman bracket of the Plat traced 
braidings. 




FIG. 4: The Markov (a) and Plat (b) tracing of a braid. In (b) the equivalence of the Plat and the Markov traces is depicted, which 
is established with the help of the Temperley-Lieb algebra element <I> 1 42] that has a simple diagrammatic representation. 

In l42ll it has been shown that |oc)(oc| = <t , i<f , 3...<f , „_i/fl!" / ' 2 . Intuitively, this means that the density matrix of the 
initial anyonic state, given by the pair creation of particles from the vacuum, corresponds to the product of Temperley- 
Lieb algebra elements <I>, as depicted in Fig. [4] Thus one has 

tr(fi>)(a|^, + ) = _i_ r ((2^V 3) Ptat )j 

where Bl,B^ have been Plat traced and then the corresponding Kauffman bracket has been evaluated. The Kauffman 
bracket 14311 of a link (L)(A) is a Laurent polynomial in the complex variable A and can be computed by smoothing 
each crossing according to the Skein relation shown in Fig. [5] Up to a phase factor they are equal to the Jones 
polynomials B44I1 . topological invariants that facilitate the distinction between inequivalent links. Kauffman brackets 
have the following property, (LO) = d(L) for L being a general non-empty link and O being the trivial link. Also 
(O) = 1. Hence (OO...O) = d'"^ 1 where the bracket is over m trivial links that correspond to n = 2m anyons. In 
particular, when B$ and By are trivial (corresponding to trivial statistics or involving trivial anyonic braidings) then 
tr(5t|a) (a\B' s , T ) = 1, i.e. we obtain the usual quantum walk, with no anyonic effect appearing. 
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0— ► -(A- 2 + A^)^d 

FIG. 5: Method to recursively calculate the Kauffman bracket (L) (A) of an unoriented link L. All crossing in a projection of the link 
can be smoothed using the Skein relations and components which are unknots contribute a multiplicative factor of d. Normalization 
is chosen so that (O) = 1. 



A. The SU(2) level k anyonic model 



While the main structures presented here are applicable to general non-Abelian models we present, for concreteness, 
analytic calculations for the SU(2)a: topological models. These have A = i e m l 2 ( k + 2 ) so their quantum dimension is 
d = 2cos(jq3) where k is an integer greater or equal to 2. All the braiding and fusion matrices of these models are 
well known 113 611 . In the following we explicitly calculate the probability distribution in the position of the walker for 
t = 1, 2, 3, 4 and 10 number of walker steps. 

Let us consider the t — 1 case. For that walk we use two pairs of anyons canonically ordered on a line and the walker 
is the third one from the left. Then we have 

= Oldiag = |2><2|trZ7^ ti-9oo+ |4)(4|tr« 1 1 1 tr^ 1 1 1 . 

ci. ci' 

But 

«ob = «ii = md fr^bo = tr 9n = ;j (OO) = 1, 

so 

£p flV (f = l)| diag = II(00)(|2)(2| + |4)(4|) = I(|2)(2| + |4)(4|). 

a,a f 

Hence, after one step the walker has equal probability (1/2) to be at the two neighboring sites for all types of anyons 
and no anyonic effects appear. 

Let us consider the t = 2 case. Here we use three pairs of anyons, where the walker is initially located at position 
so = 3. One can evaluate 

2 1 

= 4 

for all relevant pairs aa! , i.e. (0,0) (0,0), (1, 1)(1, 1), (0, 1)(0, 1) and (1,0)(1,0) while it is zero otherwise. For the 
trace of braidings we have 

^(0,0X0,0) = ^{(bibi)%b 2 ) - ^(OOO) = 1 

and similarly 

tr ^(i,i)(i,i) =tr ^"(o,i)(o,i) = tr ^(lo)(i,o) = !• 
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Hence, we finally have 



= 2)|diag = + |5)<5|) + i|3)<3|. 



Again the anyonic character is not apparent at t = 2 either. 



Let as consider the probability distribution for t = 3 where we use 5 pairs of anyons. The paths aa' that start at the 
common initial position sq = 5 and end at the same final position having at the same time the same final coin state 
are given by (000)(000), (001)(001), (010)(010), (100)(100), (011)(011), (101)(101), (110)(110), (111), (111), 
(010)(100), (100)(010), (011)(101),and (101)(011). One can see that for all of the paths we havetrwl, = 1/8 apart 
from the last two (01 1 ) (101 ), and (101)(01 1) for which tr«L, = -1/8. Also for all of the paths we have tr^l, = 1 
apart from the last four for which it is 

at -2 at i 1 2 COS ( ) COS ( t^) 

tr(B 010 ) V)(«l« too)) +fr(4oo) l«X«l B (oio)) = ~A{L) + (L)*) = ^ ^ k+2 



and 



at a it i 1 2cOs(t^ct)cOs(t^t) 

tr(B 3 (on ; |a)(a|4 01) ) +tr(^ 101) t |a)(a|^ 011) ) = -«L) + (L) ) = ■ " " 



where L is the non-trivial link given in Fig. [(J a). Hence, the position distribution is of the form 



1 /3 1 cos ( -pf^) cos (t^pt) , 

L Pa,a> (t = 3)|diag = i |2) <2| + f I + i U + 2 J s ^ ) |4) (4| - 



4 cos(^) 



3 lcos(^)cos(-^h 
' " ,;! - k+1 , |6)(6| + ^|8)(8|. 



8 4 cos(^) 



1 



The trace of this reduced density matrix is 1 as expected. Nevertheless, it acquires nontrivial contributions from the 
anyonic nature of the walker and its environment. 

For the t = 4 case we have 



E 9r,3 (t = 4) | diag = — 1 1) (1 1 + - ( -3 cos 



2jc 

k + 2 



3 cos 



4% 

k + 2 



5 1 13><31- 



1 / / 271 \ / 4tc A / 6tc 

— 3 cos — cos — 3 cos 

32 V U + 27 \k + 2j U + 2 



5 sec 2 



k + 2 



|5)(5| 



_L 2 „ s (*L - cos +1 W ^ „,(7 l + ± lm 



where the contributing links are given in Fig. [6ja-d). 




(b) 



(c) 



(d) 



JO ^ ^ 



FIG. 6: The links that contribute to the t = 3 (a) and t = 4 (a-d) anyonic quantum walks. 



We have obtained analytic results for the distribution of the walker for a t = 10 step walk. This distribution is 
depicted in Fig. Qa) as a function of the position of the walker and the level k of the SU(2)t theory. We notice that 
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for large k the distribution, P(s,t), approaches an asymptotic value. This is true for any t we checked. In order to 
determine the behavior of the walker distribution we compare it to the classical and quantum walks. For that we 
introduce the distance between the anyonic distribution P(s, 10) and the corresponding standard quantum or classical 
ones Pq,c(s, 10) as d q , c = \/Y,s \P( S , 10) — P q ,c{s, 10) | 2 . From Fig. |7Jb) we see that d q goes to zero for large k. For 
example, for ^ > 70 it is d q < 0.01. Hence, for large k the anyonic walk distributions are very close to ones obtained 
from the corresponding quantum walks. This is in agreement with the intuition that for k — ► °° the anyons behave 
as spin- 1/2 fermions (0 — 7C of the Abelian case), which reproduces the standard quantum walk distribution. On the 
other hand the distance to the classical random walk, d c , is closer to zero for small k, while in this region d„ becomes 
maximal. From Fig. |7Jb) we see that the minimum of d c = 0.126 is obtained for k = 6. Clearly, walks with much 
larger number of steps is needed in order to conclusively determine the t asymptotic behavior of this anyonic model. 




FIG. 7: (a) Distribution of the quantum walker, P(s,t = 10), for a t = 10 step walk as a function of position s and the type of 
anyonic model indexed by k. The origin has been shifted so that the starting point is i'o = 0. (b) The distance d q and d c between 
the distribution P(s,t = 10) and the corresponding quantum Pq(s, 10) and classical P c (s, 10), respectively. Clearly, for small k the 
dispersion is closer to the classical one and for large k it is closer to the quantum one. 



It is instructive to consider the decomposition of the anyonic Hilbert space in terms of qubits. As a specific example 
consider the case of SU(2)2 anyons with labels l,fJ,\|/ representing particles with spin 0, A, 1 respectively. The non- 
trivial fusion rules are 

oxo = l+f, ax\|/ = a, \|/x\|/=l, 
and the relevant non-trivial F and R matrices are 



1 



expressed in the basis {l,y}. For the total charge zero superselection sector, the fusion rules require that the only 
physical fusion basis states are those with a?/t = for k = [1, S — 1] flZ. This then gives a fusion space isomorphic 
to m = I — 1 qubits. The generators of the braid group assume a particularly simple representation in ^fusion (1) — 

(C 2 )® m 

Bi = R®y =2 h, B 2 =B®y =2 i 2 , B 3 =A®y =3 i 2 

B 2 k = ® k jZ\h®B®y =k+l l 2l Bik+i = ®yl}l2<S>A<g>^ +2 l2; \<k<m 
B„-3 = ®p! 2 l2®A, B„-2 = ®7= 1 1 l2 8)B, Bn-i^p 1 ^®*, (17) 



where 



R ~ ( ) ' B 




B. The D(G) quantum double model 
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So far we have described quantum walks with non-Abelian anyons corresponding to irreducible representations 
(irreps) of the quantum group SU (2)%. Yet there are other physical theories with anyonic excitations that are described 
by finite dimensional irreps of quasi-triangular Hopf algebras. As with quantum groups it is possible to find represen- 
tations of the braid group induced from these algebras which then allow us to compute the entanglement of the fusion 
degrees of freedom with the walker's spatial degree of freedom. First we write down a general link polynomial L(z,z) 
which is a Laurent polynomial in a complex parameter z and its conjugate z for oriented links. We demand that L(z,z) 
is an invariant under ambient isotopy. Writing the n component link L = (g) Markov a s the Markov trace over a braid 
word B = Y\ k b'p k k G <B n , with p k G {1,2, . . . ,n - 1}, m,- G Z , we have ll46ll47ll 



L{z,z) = {zz)- n l\-) <p(fi), 



where e(B) — Y,k m k (this is the same as the writhe when we orient all components in the same way). Here the quantity 
<p(Z?) is the Markov trace which satisfies the following properties 

(I) (p(fifi') = (p(fi'fi) MB,B' G <B„ 

(II) <p(Bb n )= Z<S?(B) V5G®„-iC«„, 
<p(Bb- 1 ) = z(p(B) MB G s„-i C «„, 

where z = ty(b n ), z = ty(b~ l ). As a consequence of these properties, the link invariant satisfies 

L{BB') Mmk0V = L(B') Mmkm \JB,B' G <B„, 
L(Bb^) Ma[kov = L(fi) Markov VBe ViCS„. 

Using the additivity of the writhe number it can been shown that for the link invariant L„ constructed from (p„ we have 
that, L„(0i#02) = L,(9i)L„_,(02), where L,(.) and L„-j(.) are respectively constructed from cp,(.) and cp„_,(.). 

There is a procedure to construct the braid group representation of quasi-triangular Hopf algebras but we focus here 
on a large class given by the quantum double of a finite group D(G). The irreps of D(G) are labelled by the pair 
(k,a), where k labels a conjugacy class Ck of G and a labels an irrep, with dimension |oc|, of the centralizer of a 
representative element gk G Q. The quantum dimension of these irreps (our anyons) is d[k,a] — \ Ck\\&\- From the 
representation theory of quantum doubles it is shown [48, 49] that for the invariant above 

L(z,z)=d[k,a} n {gk)a e{B) <?(B), 5£8„ (18) 

with 

m(m _ tr ( g ) (8k) a (8k l )a ,±u _ Xq(gf ) 

Vl ' {d[k,a})"+ 1 ' d[k,a]' d[k,a]' KSk >a |a| ' 
and Xa is the character in the representation a. 



Consider as an example the symmetric group G = Sn, N > 5, and C2 C Sn the conjugacy class of transpositions of 
Sn , of order \Cz\ = jN(N — 1 ) with chosen representative g2 — (N—1N), and let Z2 = Sn-2 x S2 be the centralizer 
subgroup of g2 where ^2 = {/, g2]- Let a = 1 labelling the trivial irrep of Z2, from which the (k, a) = (2, 1) irrep 
of D(Sn) with dimension d[2, 1] = lN(N — 1) is induced (see [48]). For this irrep (#£=2)01=1 = 1 an d •HU becomes 



L(B) 



Markov 



:</[2,l]"q>(2?). 



Consider the braid group on n + 1 strands, ®„, and a canonical word B = (b^)" 11 (bj 7 ) mi ... {bi„) m " G ®„. If 
(ii,i*2) ■■■,in) is an arbitrary permutation of (1,2, ...,n) and m, £ Z, the link polynomial for the trace closure of these 
words reads 



L(B) 



Markov 



= ri{l + f(iV-2) 



l+2cos(^— ) 



+ l(N-2)(N-3)[l + (-ir] 



(19) 
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By virtue of ( fTSI l. and by taking for coin initial state \\\t), we have for the diagonal elements of the walker's density 
matrix the following expression 



t j t k 

LPa,«'(0|diag = E E L E 5 *,V/S w ,Vi «^ 

2,2' ;=l^=0A-=lv { .=0 



( P((n( § ^.^-l +V,i« r -l+l)^+^+^-l-r)(n( 5 Vr-lAr+ 8 Vr-l,Vr-l)^+Vr+Vr-l-r) + ) 



r=l r=l 

(20) 



n+1 „ \/n+l „ 
— — +2(i,-t )( — — +2/i f -/ 



where « ftV = tr (n' r= o(^,+ 1-^,^)1 ¥)(¥ I (lT=o(Av +1 -v,.f/)) J. Here we take = |0)(0| and p walker = 

|t~7 (^tH as tne ™^ com anc l wa lker density matrices, respectively. As we now employed the Markov trace 
rather than the Plat one the obtained probability distribution corresponds to the modified quantum walker algorithm 
where we start with n anyonic vacuum pairs, braid all right members of each pair uniformly to the right and then per- 
form the quantum walk just on the n strands corresponding to these right members (see Fig.Ua)). Also the reshuffling 

matrix is taken to be U = -4= ( \ \ ] . It is worth noting that the same results reported below for the U reshuffling 



</2\i 1 f 

matrix are valid if we had employed the Hadamard matrix H, as a straightforward calculation shows. 

Next we briefly present the expression for the diagonal part of the walker's density matrix for the first time steps. 
For step t = 1 , we obtain 

I P2,2'(* = Old* = \v{B\B 2 ) |2) (2| + ~q>(B^ 3 ) |4) <4| = ~ |2) (2| + ~ |4) <4| . (21) 

a, a' 

For step t — 2 we have 

X>2,3'(' = 2)|diag = 1 -{(?((B 1 B 2 ^B 1 B 2 )\l) (l|+(p(( J B 4 B3) t B4 J B 3 )|5) (5| 

a,a' 

+ [c P ((B 2 ) t B2) + q>((fi 3 ) + fi 3 )] |3> <3| } 
= I|l)(l| + i|3)(3| + l|5)<5|. (22) 

We notice that the diagonal part of the density matrix given above for the first two steps is common (i.e. independent 
of N), to all cases of the symmetric group used i.e. Sn, N >5. This changes for subsequent steps as we see below. 

For step t = 3 we calculate the expression 

Y,Pa,a'( f = 3)| diag = ^{|2> (2| + (3 + 9 ((Bl)^B 4 ) +9((B|B 4 ) t B|)) |4) (4| 

a, a' 

+ (3 - cp((4B 5 )^) - <?((BlMB 5 )) |6) <6| + |8) <8|} 
= I{|2) (2| + (3 +<v(BfB 2 3 ) + <p(B^ 3 " 2 )) I 4 ) ( 4 I 

+ (3 - ^{B\B 6 2 ) - «?{B- 5 2 Bl)) |6) (6| + |8> (8|}, (23) 

or in terms of link polynomials 
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£piMf (* = 3)|dia g = i{ |2) (2| + (3 + J[2, 1]- 9 L(B 4 2 Z?2) +£ /[2, l]- 9 L(fl 2 Z? 3 - 2 )) |4) (4| 

+ (3 - d[2, 1]- 9 L(8^ 6 2 ) - d[2, 1}- 9 L(BJ 2 B 2 6 )) |6) (6| + |8> (8| } 
= i{|2)(2| + (3+2,/[2,l]- 2 [l + i(^-2)(^-3)] 2 )|4)(4| 

+ (3-2,/[2 ; l]- 2 [l + i(^-2)(^-3)] 2 )|6)(6| + |8)<8|}, (24) 
where (fT8l ) has been used. By means of ( fl9] l we obtain 

L(B 2 B 6 2 )=L(BfB 2 )=L(B 4 2 B 2 )=L(B 2 Bf) = [l + 2)(tf- 3)] 2 d[2, l] 7 

For the choice N = 5, the dimension is af[2, 1] — 10, and this yields 

P(' = 3)1^ = \{ |2) (2| + (3 + 1) |4> (4| + (3 - 1) |6> (6| + |8) (8| }. (25) 

Similarly for the f = 4 step we obtain 

LP^'C = 4)| diag = l|l)(l| + {i + ^[l + i(^V-2)(^V-3)] 2 ,/[2,l]- 2 }|3)(3| 

a : a f 

+ {^-^[l + ^(^-2)(/V-3)] 2 £ /[2,l]- 2 }|5) <5| 

+ {j6 W [1 + ^(^- 2 )(^- 3 )] 2fl, [2' i]" 2 } I 7 ) (71 + ^ 19) <9| , (26) 
which for the N — 5 and d[2, 1] = 10 give the diagonal part 

£ Pa^(* = 4)|dia g = ^ |1) <1| + ^ |3) <3| + |5) <5| + ^ |7) <7| + 1 19) (9| . (27) 



Before closing this section some comments are in order. Comparison of the diagonal parts of walkers's density 
matrix obtained above for the two forms of modelling the braid group representations, i.e. quantum group SU (2)t 
and symmetric group respectively, shows that the walker occupies the same sites at both cases albeit with different 
occupation probabilities. This difference in the distribution of occupation probabilities is expected to be a general 
feature, since by construction these probabilities depend both on the algebraic structure (finite group, quantum group 
etc.), employed to provide a matrix representation of the braid group, via the associated link polynomials. It is worth 
mentioning that in the limits k — > °o and jV — > °°, both distributions converge to the standard quantum walk as seen in 
Fig.! 



V. OUTLOOK 



We have described here a simple quantum walk on a line segment with n sites. This is a path graph with vertices 
given by the locations of n anyons and n — 1 edges connecting the vertices that correspond to generators of the braid 
group B„. One could imagine anyonic walks on more general graphs by embedding the graph onto a cellular surface of 
genus g. Such a mapping is not unique and the study of such embeddings is the focus of topological graph theory [50]. 
One could then perform walks with an anyon on this higher genus surface using the representation theory for braid 
groups on arbitrary manifolds 115 ill . 
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FIG. 8: Combined surface plot for the distribution of a t = 4 walk in the modified protocol that corresponds to the Markov trace. 
The blue surface is the distribution for spin- 1/2 anyons in the quantum group SU{2)i ( and the green surface for quantum double 
anyons in the (2, 1) irrep of D(Sn). 



There are numerous questions on anyonic quantum walks that have not been addressed yet. The Abelian anyonic 
quantum walk is easily tractable, both theoretically and numerically, giving rise to a distribution similar to the standard 
quantum walk. As we have seen one could determine the behavior of non- Abelian quantum walks for a small number 
of steps. It is a very hard problem to probe any non-Abelian walk for large anyonic chains due to the exponential 
increase in their Hilbert space dimension. Each additional step of the walker exponentially increases the number of 
different paths that one needs to consider. On the top of that one needs to evaluate the Kauffman brackets of the 
corresponding paths in order to determine the final position distribution of the walker. For different types of anyonic 
models (like SU(2)2) the evaluation of the Jones polynomials is computationally polynomially demanding lf52l l53ll . 
though still one has to consider the exponential increase in the number of paths. For other models (like SU(2)i for 
k > 2 and k ^ 4) the Jones polynomial evaluation is believed to be, in general, exponentially hard to compute I53l 15411 . 

Ideally, one is interested in deriving the asymptotic behavior of non-Abelian anyonic quantum walks. The examples 
we considered here reveal parts of the structure and possibly indicate the asymptotic behavior of certain models. For 
example, we explicitly considered the SU(2)^ anyonic model for up to t = 10 number of steps. This gave the following 
interesting characteristics. For small k the walk appeared to have a distribution close to the classical random walk. For 
large k the distribution appeared to be similar to a quantum. Even though small t distributions can hardly be conclusive 
about the large t behavior they provide the exciting possibility of having an anyonic model with intermediate values of 
k where the variance is in between the classical and the quantum. That would constitute a new paradigm of diffusive 
behavior. A variety of non-Abelian models needs to be systematically considered theoretically and numerically in 
order to conclusively understand the asymptotic behavior of non-Abelian walks. 
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